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ON THE STABILITY OF SOLUTIONS OF SEMILINEAR ELLIPTIC EQUATIONS 
WITH ROBIN BOUNDARY CONDITIONS ON RIEMANNIAN MANIFOLDS 


C. BANDLeQ, P. MASTROLIi^, D. D. MONTICELlI^ AND F. PUNZcQ 

Abstract. We investigate existence and nonexistence of stationary stable nonconstant solutions, i.e. 
patterns, of semilinear parabolic problems in bounded domains of Riemannian manifolds satisfying 
Robin boundary conditions. These problems arise in several models in applications, in particular in 
Mathematical Biology. We point out the role both of the nonlinearity and of geometric objects such as 
the Ricci curvature of the manifold, the second fundamental form of the boundary of the domain and 
its mean curvature. Special attention is devoted to surfaces of revolution and to spherically symmetric 
manifolds, where we prove refined results. 


1. Introduction 


In this paper we study stability and instability of solutions of 


( 1 . 1 ) 


du 

Au + fiu) = 0 in n, ——h art = 0 on dVl. 
ov 


Here H is a smooth domain in a Riemannian manifold (M, g)^ A is the Laplace-Beltrami operator on M, 
V is the outer normal to with respect to M and a £ R is an arbitrary fixed number. A solution of 
problem CH) may be regarded as a stationary solution of the parabolic problem 

du 


( 1 . 2 ) 


dt 


= Au + f{u) in H X (0, T) 


du 

-—\- au = \) 
dv 


in dVt X (0, T) 


u = uq in H X {0} . 

Stable solutions are used in mathematical models of pattern formation. They are often called patterns 
and have attracted much attention in the literature (see e.g. [5], [3], [1], [S], [H], [S], [S], [IS], [S], 
[22) . [25)1. Applications of problem (HI) to mathematical biology are found e.g. in ig, HD, im. A 
biochemical process on surfaces of revolutions is described and analyzed in m- In most papers it is 
assumed that the boundary is impermeable (that is a = 0). However it is reasonable to consider also 
the case of a flux that is proportional to the solution (see e.g. |g, [ig, jUj). This motivates our choice 
of Robin boundary conditions with a £ R. 

We recall (see e.g. (ig) that a solution of problem (11.11) is asymptotically stable if the smallest 
eigenvalue Ai of the linearized problem 

Afj) + f'{u)(f> + Xcj) = 0 in O 

(1-3) ! 

90 , . 

——\- ag) = 0 mail 

ov 
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is positive, it is unstable if Ai is negative and neutrally stable if Ai = 0. Asymptotically stable solutions 
U of problem (EU have the property that they attract for large time all solutions of E3 which initially 
are sufficiently close to U (see e.g. [TO]) while the latter are repelled from the unstable ones. If Ai = 0 
both situations may occur. 

The sign of Ai is also crucial in the calculus of variations. Indeed, let {xi}"^ be a system of local 
coordinates and let gij be the corresponding metric tensor of M. Its inverse will be denoted by gb 
Furthermore we have g = g''^Uxi<l)xj and |VmP = (;(Vu, Vu). Here and in the sequel we 

shall use the Einstein summation convention. The solutions of (HH) are related to the following energy 
functional 

£(u,fl) = f iVupd/r + a/ dS — 2 f F{u) dg, 

J Q J dQ J O 

where dfi is the Riemannian volume element, dS is the surface element of 3H, and F'(u) = f{u). If u is 
a solution of (HE then the Frechet derivative of £ vanishes at u, more precisely 

—£{u,£l)= I g {Vu,\/cj)) dfi + a (p uipdS— j f{u)(j)dg = Q for all 
2 Jn JdQ Jn 

The second derivative of f at u is 

^£{u,£l) = ( \V(j)\^ dg + a cj)^ dS — ( f'{u)(j)^dg. 

^ J Q J 50 J O 

If Ai is positive, then £{u, £1) is a local minimum whereas if Ai is negative u is a saddle point. 

The study of stability of the solutions of (HU has a long history. First results were obtained by 
Hudjaev m and Keller and Cohen m for problems in K.", while the stability and the uniqueness of 
solutions to (HH) in the case of weighted concave nonlinearities and positive solutions were studied in 

E). 

Casten and Holland |1], and also Matano [TO], observed that for problems with Neumann boundary 
conditions (a = 0) in a convex domain in R", all nonconstant solutions are unstable. This result was 
generalized to problems on manifolds by Jimbo [TO] and by Handle, Punzo and Tesei [2] (see also [2T]b 
Jimbo proved that all non-stationary solutions are unstable if the Ricci curvature of M and the second 
fundamental form of the boundary are positive. The aim of this paper is to study the stability of solutions 
to problems with Robin boundary conditions in bounded domains, both in R'” and, more in general, in 
Riemannian manifolds. In this case the conditions on the boundary, which imply instability, depend also 
on the nonlinearity /. 

To give an idea of our results, let 17 be a domain in R"*, let i = 1, ..m — 1, be the principal 
curvatures of d£l and denote by FI = — (m — 1)“^ mean curvature of d£l. By Theorem 14.41 

we have that ifa + (m — l)iF + < 0 on dfl and if an additional assumption involving the second 

fundamental form of dfl and a is satisfied, then every non trivial solution is unstable. Note that this 
condition in case m = 2 reduces to a — IF > 0. For non positive a in particular it implies that dQ is 
convex. 

As a counterpart we can derive from this type of considerations estimates for stable solutions. It 
should be pointed out that no assumption on the sign of a or on the solution is made. In the case of 
Riemannian manifolds a similar result holds under the additional assumption that the Ricci curvature is 
nonnegative. However for surfaces of revolution or for problems on spherically symmetric manifolds we 
can allow the Ricci curvature to be negative provided it satisfies a suitable bound from below, see Sections 
|5|and|6l Furthermore we construct by means of arguments developed in [2] and in [2^ a counterexample 
which shows that this bound is sharp . 

The discussion of nonexistence of stable solutions is based on the variational characterization of 
Ai, on the well-known Bochner-Weitzenbock formula and on a, to our knowledge, new result on the 
decomposition of the normal derivative of |Vr(;(a;)p on dfl where w satisfies di,w = —aw on dVl, for some 
a G R, see Theorem 13.41 Similar formulas are known in the literature for special cases and have been 
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extensively used by L.E. Payne and his collaborators to obtain estimates for the solutions of boundary 
value problems see e.g. [53]. 

The particular case of Neumann boundary conditions a = 0 was studied in [T3] for 17 C R™ and in 
[2] for 17 C M. The case of a general a G R involves a deeper rather technical analysis based on the 
method of moving frames, see Section |3| 

The paper is organized as follows. At first in Section [51 we deal with the case of domains of in 
order to illustrate the flavor of the technique. Section|3]is an introduction to the geometrical notions and 
tools needed in the sequel. Section [4] is concerned with the general result on instability on Riemannian 
manifolds, while Section |S] is devoted to the case of surfaces of revolution where sharper results are 
obtained. The last section contains some applications to specific manifolds. 


2. Domains in the plane 


In order to get a better insight and because the arguments are elementary we treat first the case where 
17 is a domain in the plane. 

Suppose that the boundary 917 is represented by the curve s i-A x{s) := (xi(s),a; 2 (s)), with s G [0,/], 
where s is the arc-length. We assume that x{s) is positively oriented, and that the curve is sufficiently 
smooth, so that the differential equation (HH) holds up to the boundary. The outer normal to 917 will 
be denoted by v. In a neighborhood of the boundary a point cc G 17 is given by 


( 2 . 1 ) 


x{p, s) = x{s) — pv{s). 


Here (p, s) are called normal coordinates. By the Frenet formula v = nx where k is the curvature of 917. 
The metric g can be written as p = (1 — pK)^ds'^ -I- we thus have 


9ij — 


(1 — pk)^ 0 


0 


1 


= 


(1 — pn) ^ 0 


0 


1 


0 ^ = det(gij) = (1 -pk)^ 


where pij are the components of g and are the components of its inverse. Moreover, the volume element 
dx is dx = {1 — pK)dsdp, while for a sufficiently regular function u the gradient and the Laplacian (see 
also the next section) are, respectively, given by 

»,2 


( 2 . 2 ) 


|Vu|2 = 
Au = 


(l-pit)2 

1 / u. 


1 


-{{1- pK)Up)p. 


1 \ 1 / ' 1 r—/'^pp- 

\ — pn \ \ — pn J ^ \ — pn 

Consider now the original problem (HU and the corresponding eigenvalue problem HU. The smallest 
eigenvalue is characterized by the Rayleigh principle 


(2.3) 


Ai = min 


/n dx - 4 f{u)v^ dx + a§. 


an 


ds 


In dx 


Following Casten and Holland |2], we choose v = Uxi as a test function in (I531) . Then by the Gauss 
theorem 

Ai / u^. dx< dx + a (p u^. ds — f'{u)u1. dx 

JQ ' Jn Jan ' Jn 

= — I UxiAuxi dx + - (p ds + a (p ds — j f'{u)u^. dx. 

Jn 2 Jqq ' Jqq ' Jq 

Replacing Auxi by —f'{u)uxi and summing over i we get 

(2.4) Ai I |Vup dx < i ® d^\Vu\'^ ds + a (p |Vup ds. 

Jn 2 Jqq Jqq 
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Next we compute the first integral on the right-hand side of (12.41) . From (12.21) we conclude that on dfl 
(that is p = 0 in ( 12 . 11) 1 

= -^dp\Vuf = -UsUsp - Kul - UpUpp. 

Keeping in mind the boundary condition Up = au and from ( 1 ^ and (HH) the relation 

Upp — KUp + Uss + f{u) = 0 on 


we obtain 
(2.5) 

Since 


= —{a + k)uI — Ka^u^ + auuss + auf{u). 


a ® uuss dS = —a ® uZ ds 


an 


an 


we finally get 

(2.6) Ai / |Vm|^ dx< (f 

Jq Jd 


2 2 
a u 


an 


a — K-\- 


fiu) 


— (a -|- k)Ms ds. 


Theorem 2.1. Assume a ^ 0. If on dil the conditions 


(CO) / 

Ja 

(Cl) 

are satisfied then Ai < 0 and u is unstable. 


2 2 
a u 


a — K,-\- 


fiu) 


an 

O. ^ 0 - 


< 0 


Observe that 

(C2) a — Kmin -|-< 0 

au 

implies (CO). 

Remark 2.2. If dil is a level line of u, then Ug = 0 on 90 and hence Theorem 12.11 remains true only 
under the condition (CO). 


If 90 consists of several closed curves Fj, i = then condition (C2) can be replaced by a — 

minr^ k + maxr^ < 0 for i = 1,..., k. 


Example 2.3. Let O be an annulus with the radii tq < R and u = u{r) be radial. Then u is unstable if 


roa -I-1 -I- 


rofiu{ro)) 


^ifo) 


M^(fo) 


Ra — 1 


RfHR)) 


i(i?) J 


u^iR) < 0 . 


Note that inequality (EH) contains also information on stable solutions. In fact for a stable solution 
Ai is positive and thus 


0 < 


an 


2 2 
a u 


a — K + 


/(w) 


au 


— {a + k)i 


ds. 


If (Cl) holds then the expression a — k+ must be positive somewhere. For instance if u is the first 
eigenfunction of Au -|- Am = 0 in O, ^ om = 0 on 90 and a > 0 then 


Ai ^ CH . 


3. Some useful tools from Riemannian geometry 

In this section we collect some notions and results from Riemannian Geometry following |18j and [T] . 
Moreover we prove a decomposition theorem (see Theorem 13.41 below) for the normal derivative of the 
squared norm of the gradient of an arbitrary smooth function satisfying Robin boundary conditions. 
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3.1. Basics on the method of moving frames. Let {M,g) be a Riemannian manifold of dimension 
m with metric g. Let p £ M and let {U, ip) be a local chart such that p £ U. Denote by ..., x™, 
m = dim M the coordinate functions on U. Then, at any q £ U we have 

(3.1) g = gij dx''® dx^ 

where dx* denotes the differential of the function x* and gij are the (local) components of the metric 
defined by gij = g (^r, gfr)- In equation (13.11) and throughout this section we adopt the Einstein sum¬ 
mation convention over repeated indices. Applying in q the Gram-Schmidt orthonormalization process 
we can find linear combinations of the 1-forms dx* which we will call d® for i = 1,... ,m. Then (EH) 
takes the form 

(3.2) g = 5,je^®e\ 

where Sij is the Kronecker symbol. Since, as q varies in U, the previous process gives rise to coefficients 
that are C°° functions of q, the set of 1-forms {d®} defines an orthonormal system on U for the metric 
g, i.e. a (local) orthonormal coframe. It is usual to write 

m 

instead of (EH- We also define the (local) dual orthonormal frame {ci}, for z = 1,... ,m, as the set of 
vector fields on U satisfying 

(3.3) 0He.)=Si, 
where Sf is the Kronecker symbol. We have the following 

Proposition 3.1. Let {d®} be a local orthonormal coframe defined on the open set U C M; then on U 
there exist unique 1-forms {d® }, for i, j = 1... ,m, such that 


(3.4) 

dd® = -d® A d^ 

and 


(3.5) 

d® + d^ = 0 


The forms d® are called the Levi-Civita connections forms associated to the orthonormal coframe {d®}, 
while equation (El is called the first structure equation. 

The curvature forms {©}} are associated to the orthonormal coframe {d®} through the second structure 
equation 

dd® = + Q]. 

Because of (13.51) it follows immediately that 

(3.6) 0® + 0^ = 0. 

Using the basis {d® A d-^ }, for 1 < i < j < m, of the space of skew-symmetric 2-forms h.^{U) on the open 
set U, we may write 

(3.7) 0® = ii?®fc,d^ A d* 
for some coefficients i?®^.^ G C°°{U) satisfying 

(3.8) R®,,+i?®,, = 0. 

These are the coefficients of the (1, 3)-version of the Riemann curvature tensor which we denote by R. 
More precisely, in this local orthonormal frame we have 

R]kt = &]iek,et) = (dd® -f dfc Ad^)(e/c,et) = g{R{ek,et)ej,ei), 
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SO that its components are 

R = ® 9^ ® e*. 

Note that (13.61) implies 

(3-9) ^}a+^L=0- 

The (0,4)-version of R is defined by Riem(Jf, T, Z, bh) = g{R{Z,W)Y, X), so that its local coefficients 
Rijkt satisfy 

Rijkt = Riem(ei, ej, Cfc, e*) = g[R{ek, et)ej,ei) = R]kt 
and thus in the local orthonormal frame 


(3.10) Riem = Rijkt9^ ® 9^ ® 9^^ ® 9*. 

For further details, we refer to [T]. The Ricci tensor is the symmetric (0,2)-tensor obtained from 
(13.1011 by tracing either with respect to i and k or, equivalently, with respect to j and t. Thus 

Ric = Rij9^ ® 9^ 

with 

Rij — Ritjt — Rtitj- 

Now let u S for the differential of u, du, we can write 

(3.11) du = Ui9^ 

for some smooth coefficients Ui\ the Hessian of u is then defined as the (0, 2) tensor field Hess(M) = Vdu 
of components Uij given by 

(3.12) Uij9^ = dui - Uk9i, 
that is 

(3.13) Hess(M) = Uij9^ ® 

Here and in what follows V is the (unique) Levi-Civita connection associated to the metric g. It is easy 
to prove that 


(3.14) — ‘^jit 

so that Hess(it) is a symmetric tensor. In global notation we have, for all smooth vector fields X, Y on 

M, 

Hess{u)iX,Y) = {Vdu){X,Y) = Y{X{u)) - {VyX){u) = X{Y{u)) - (Vxr)(u); 
it is also possible to show that, equivalently. 


(3.15) 


Hess(u)(X,r) = ^{Cvug){X,Y), 


where C'^uQ is the Lie derivative of the metric g in the direction of Xu. With respect to local coordinates 
{a;®}, i = 1,..., m, the Hessian is given by 


(Hess(u))„. = 


d'^i 


-rf 


du 


dxidxi ^^^dxk 


d ^ d 
ssj dx, ~ si- dx, ■ 


where T^^ are the Christoffel symbols, defined as usual as 

d 

dxk "“"J 

In the moving frame formalism the squared norm | Hess(u)p is given by UijUtj, while in coordinates we 
have 

|Hess(u)|" = gV' 


/ d'^u 

rk 9u\ 

/ d^u 

r« 9u\ 

\dxidxj 

1 

\dxkdxi 

''' dxs) 
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The Laplacian of u is, by definition, the trace of the Hessian, (more precisely, of the (1,1) version of 
the Hessian), that is, 


Am = Tr(Hess(M)) = uu. 


The gradient of a function m : M —R relative to the metric of M, Vm, is the vector dual to the 1-form 
du, that is 

g{yu,X) = du{X) = X{u). 


for all smooth vector fields X on M. In a local orthonormal frame we have Vm = u^Ci = UiCi, so that 
|Vm|^ = UiUi, while in local coordinates we have 


\Vu\^=g 


du du 
dxi dxj 


The divergence of a vector field V = V^Ci on M is given by the trace of W, the covariant derivative of 
H; since XV = (dH* -I- ® Ci = V-9^ 0 Ci we have 


divH = 


while in local coordinates 

dV^ ■ 

Note that the Laplacian of the function u is the divergence of its gradient, that is 

Am = div (Vm). 


In local coordinates it has the form 

The third derivatives of u are defined by 

(3.16) UijkO'" = dutj - Ukj9f - u^kd']. 

Note that taking the covariant derivative of (13.141) we have 

(3.17) Uijk — Ujik- 
The commutation rule of the last two indices is given by 

(3.18) Uijk — Uikj UtRijk — Uikj UtRtijk- 

We state here the classical Bochner-Weitzenbock formula, see e.g. [TH], which will play a crucial role 
in our investigations. 


Lemma 3.2. For all u € C^{M) we have 

(3.19) ^A|Vu|^ = |Hess(M)|VRic(VM,VM)-k 5 (VAM,VM). 


Example 3.3. Let M C R^ be a simply connected surface and let H be a domain on M. It is well- 
known that M can be mapped conformally into R^. In this case the coordinates (xi,X 2 ) are isothermal 
and the corresponding metric tensor is gij = p^Sij. The differential operators then become 

A=p“^Afl and |Vp = Vi^p, 


where A/j and Vr are the Laplacian and the gradient in R^. In this case we have 

■UiU^kip 
f ^RU 


1 A|VmP = MifeMifeP ‘^+UiUikkP ‘^+2UiU^k{p ^)fc -f i|VfiMpAp ^ 


= UikUikP +2UiUtk{p )k-Ui{p )iXrU + Ui 


V 








Here the subscripts denote differentiation with respect to Xi. Set for short / := logp. Then (p ^)k = 
Thus in two-dimensions 

UikUrkP~'^ + 2UiUik{p~^)k - 

= p"^[(Mn -f 2 mi2 + U 22 ) + 2(lii/i - U2f2){u22 “ Un) - 4(ui/ 2 -h 1 i 2 /l)ui 2 )] 

= p"^ [(nil - Uifi -I- 142 / 2 )^ -I- {U 22 + flUl - ^ 2 / 2 )^ + 2(ui 2 - ni /2 - ^ 2 / 1 )^] 

-2p-^{ul + ul)ifi + fi). 

Moreover 

2 p^ p^ 

Note that K = —p“^A/{/ is the Gaussian curvature of M. Consequently 

i A| Vup = p~^ [(nil - ui/i + 142 / 2 )^ + {u22 + fiUi - n2/2)^ + 2(ni2 - ni /2 - n2/i)^] 

^ '-V-' 

|Hess(u)p 

+ P~'^K\Vru\'^ + p~'^Ui{Au)i. 

3.2. Immersed submanifolds. Let (A, p^v) and M be respectively a Riemannian manifold and a man¬ 
ifold of dimensions n and m, with m < n. Let / : M —>■ A be an immersion and let g = f*gN be the 
metric induced on M by / where /* denotes the pullback. If gM is a given Riemannian metric on M 
and / : M ^ A is an immersion we say that / is an isometric immersion if gM = 9 = f*Q n- 

Let V C A be an open set, and let p € /“^(V). By reducing V we can assume that the connected 
component U of f~^{V) with p is an embedded submanifold in the domain of a local flat chart. 

We fix the following indices convention: 

1 < L i, fcj • ■ • < w, m -I-1 < a, /3, 7 ,... < n, 1 < a, 6 , c,... < n. 

By means of the Gram-Schmidt procedure we can construct an orthonormal frame {Ea} in a neigh¬ 
borhood of f(U) such that {Ei} is a basis for Tf{U). We call this frame a Darboux frame along /, and we 
write {ci} for the basis of the tangent space at U such that /,ei = Ei (where /*ei is the pushforward of Ci 
by the map /). The dual {0°“} of a Darboux coframe is called a Darboux coframe along f. The definition 
of a Darboux (co)frame is equivalent to say that the vectors {Ei} (locally) span f^TM, the image of TM 
through / in TA, while the vectors {Ea} are orthogonal to f^TM and span in fact the normal bundle 
TM-^, that is the set of (local) vector fields in A that are orthogonal to f^TM. A consequence of the 
choice of a Darboux frame is that 

(3.20) f*9°‘ = 0, 

where f*0°' is the pullback of 0“ by the map /. Indeed, for every i, {f*0°‘){ei) = 0‘^{f^ei) = 9°‘{Ei) = 0. 

Let now {9^} be the Levi-Civita connection forms of A relative to {9°'}. Pulling-back on M the first 
structure equation of A, and using the properties of the pullback we have 

/*(dr) = d(rr) = -r{9^ A 9 ^) = -{r9^) a {r9^). 

By ((3.2011 we obtain in particular that 

(3.21) d{r9^)=-{r9})A{r9^); 
moreover we obviously have 

r{o}) + r{9{)=o. 

Thus from the uniqueness, see Proposition 13.11 we deduce that {f*9j} are the Levi-Givita connection 
forms of M. 

Since the pullback commutes with exterior differentiation and wedge product we shall omit from now 
on the pullback. From the context the reader should be able to distinguish between forms or tensors. 
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Then equation (13.201) becomes 

(3.22) = 0 on M 

and on M we have 

(3.23) e] +ei=Q 
and 

(3.24) = 

To obtain further information we differentiate (13.221) , use (13.241) and (13.221) again to obtain 

(3.25) Q = de°‘ = -0“ A 0* - 0^ A 0/^ = -0“ A 0b 

Hence a simple computation shows that there exist locally smooth functions such that 

(3.26) 0“ = 0J 
and 

(3.27) h% = h%. 

It can be shown that the hfj^s are the coefficients of the second fundamental tensor II : TMxTM —)■ TM-^ 
of the immersion. II is a (I, 2)-tensor along / or, equivalently, a section of T*M (g) T*M (g TM-^ (where 
TM^ is considered as a subset of the pullback bundle f*TN) which in the present setting is defined by 

(3.28) II = /i“/*(g0^(gi;«. 

One can also verify that by (13.271) II is defined globally and that it is symmetric. The mean curvature 
vector field is given by its normalized trace, that is 

H = Itr(II) = 
m m 

If is a unit normal vector field the mean curvature in the direction of v is 

/r"=5(H,0^. 

If m +1 = n and both the hyper surface M and N are orientable, we can choose Darboux frames along 
/ preserving orientations, i.e. such that 0^ A • • • A 0’"+^ and 0^ A • • • A 0"* give the correct orientations. 
More precisely the vector field Em+i dual to 0™+^ on N is, when restricted to M, a global normal vector 
field on M. We shall call it v. Furthermore note that in this case in local coordinates one has 

hi] = -g{'^eii',ej) for f,j = I,...,m, 

which in global notation can be expressed as 

(3.29) II(X, Y) = -g{yxv, Y) for any X, F e TM . 

The mean curvature in the direction of v is called the mean curvature of the immersed hypersurface 
and denoted by H. Observe that, according to our sign convention, the mean curvature of the sphere 
S™ C with respect to the outer normal is —1. Note that, with respect to the notation used 

in Section [2j if we consider 917, the boundary of a regular domain 17 C K.^, we have H = —k . 

3.3. A decomposition theorem. 

Theorem 3.4. Let {M,g) be an m-dimensional Riemannian manifold and let ft C M be a compact, 
orientable domain with boundary dfl. Let II and H denote respectively the second fundamental tensor 
and the mean curvature of the embedding dfl ^ ft in the direction of the outward unit normal vector 
field V. Let w G C^ifl); ifw satisfies 

d'UJ 

(3.30) = g (yw,v) =—aw on dfl 
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\ d ~ ~ ~ 

- —|Vu;|^ =Il(Vw, Vw) — a| Vwp — aw Hess (w)(:/, i^) on i90, 


for some a S K, then 
(3.31) 

where Vw = Vw — g (Vw, v) v is the tangential gradient with respect to dil. 

Proof. Let {ca} = {ei,..., em_i, e™ = be a Darboux frame along DM ^ M. Set 

iLy = g{ll{ei,ej),i'), 

so that 


H = 


~Hkk ; 


m — 1 

where fc<m — 1. By definition of the covariant derivative we have 


thus 


dWjfi — 

WmBd^ = WidT + dWm- 


Pulling back the previous relation to dVL and using (13.3011 we deduce 
(3.32) Wmi = HijWj — awi on dfl, 

which implies 


WiWmi = HijWiWj — aWiWi on dfl. 


(3.33) 

We now have 

thus 

19 2 1 2 

(3.34) = 2 im + WmWmm = WiWim “ aWWmm On dfl. 

Combining (I3.33P and p.34l) we get the desired result. 


i(|Vw|^)^ = WbWbA = WiWiA + WmWmA = W^WiA “ aWWmA On dfl, 


□ 


We conclude the section by recalling a relation between the Laplace-Beltrami operator A of the 
manifold {M, g) acting on a smooth function w defined in a neighborhood of the boundary 9H and the 
Laplace-Beltrami operator A of the manifold dPl, acting on the trace of the function w on dVt. Let H 
be the mean curvature of 9H and A be the Laplace-Beltrami operator of the manifold 90 endowed with 
the metric induced by the embedding 90 ^ M. Then on has, see e.g. [15] 


(3.35) 


dvj 

Aw = Aw — (m — 1)^"^—1“ Hess(w)(i^, v). 


Example 3.5. Let O be a domain on a two-dimensional surface as in ExamDle l3.3l We consider as before 
its conformal projection onto the plane. We shall use the same notation as in Section 2 and Example 
13.31 In this case we have djdv = —p~^dp and the expression (13.341) reads as 

A|Vu |2 = -p-'^\VRu\l-p~^\S/Ru\'^{p~'^)p 

Keeping in mind that p~^{k — dp logp) = Kg is the geodesic curvature of 90 we find by the arguments 
in Section 2 

1_^ 

2 9j^' 


(3.36) 


^{-UsUsp-Kul-UpUpp)+p ^pp{ul + ul) 


= -p ^{UsUsp+UpUpp) -p ^Kgul+p ^PpU^p. 


Taking into account the boundary condition Up = au we obtain 

n2 


1 d 1^^ |2 / \ Ck-U r /I \ /I \ 1 

^ - Us[logp)s - Upp\ . 
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Remark 3.6. Similar results to Tlreorem l3.4l in the special case where M = R™ were used by L.E. Payne 
and his collaborators in their study of a priori bounds for elliptic problems. This theorem provides a 
tool to determine the points where the P-function takes its maximum. A survey of these results is found 
in [23]. 


4. Instability results on Riemannian manifolds 


Let n C (M, g) be a smooth bounded domain and let u : > R be a solution of 


(4.1) 


Am + /(u) =0 in n, 
^ + au = 0 on dfl, 


where u denotes the outer normal unit vector at Oil and f € C^. Define 

In - In dgPa cj)'^ da 


(4.2) 


Ai := inf 


In 


We note that by standard elliptic theory Ai is achieved by a function (j)i € iL^(D) which is a positive 
solution of (O. 

Note that the case a = 0 corresponds to the problem of homogeneous Neumann boundary conditions 
which has already been studied in |5]. 


We first consider the case of constant solutions to problem 611 ). It follows immediately from the 
boundary conditions that tt = 0 is the only possibility. The equation implies that /(O) = 0. Let 


(4.3) Ai 
be the smallest eigenvalue of 

(4.4) 

By (|4.2I1 we have immediately 


inf 

0GiLi(a),0^O 


J^\V<j)\^dg + aJg^ 4 ,^ da 

In dn 


-Aipi = Aiipi 


D, 


+ a^pi =0 on (ID. 


Proposition 4.1. Let a 7^ 0 and let u he a constant solution of problem (14.111 . Then u = 0 and f(0) = 0. 
Moreover, for Ai defined in 63, 

i) *//^(0) > Ai then u = 0 is unstable, 

i) *//^(0) < Ai then u = 0 is asymptotically stable. 


Remark 4.2. If we use (j) = 1 € iJ^(D) as a test function in the definition (14.3p of Ai we get the upper 
bound 


Ai < a 


\dn\ 

W’ 


where 


\n\ := f dfx, 

Jn 


|5D| := 


in Jon 

This together with Proposition ld. II implies that the solution u = 


da. 

0 is unstable if /'(O) > 


We now discuss the stability of non constant solutions to problem ill. 

Proposition 4.3. Let u € be a solution of (14.11) . with f G . Then 

(4.5) Ai / \Vuf dp. < f — a\Vu\'^ + + auf{u) + a'^{m — l)u'^L['] da 

Jn Jdn ^ ^ 


Ric (Vm, Vu) dp. 
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Proof of Proposition \4-tt\ If we introduce in the variational characterization (14.21) of Ai the test function 
|Vup we get 

(4.6) Ai f \Vuf dfi < ( \V\Vu\\^dfi— ( f {u)\Vu\^ dp + a f \Vu\'^ da. 

Jq Jq J Q J dQ 

If we apply the Bochner-Weitzenbock formula (13.191) to the solution u of (14.11) and use the divergence 
theorem we obtain 


(4.7) / |Hess('u)|^ d/r “ ^ / Ric (Vit, Vu) d/r — / g{VAu,Vu)dp 

Jn 2 Jq Jq 

■ [ —iVul^dcr— f Ric (Vu, Vu) d^ + f /'(it)| Vit|^ d/r. 
Jdn ov Jn Jq 


1 

2 Jan Jn Jn 

The first integral at the right-hand side of dMl) can be estimated by means of the inequality 


|V|Vu||^ < |Hess(u)|^ 

This result follows immediately from Schwarz’s inequality. Indeed if we use a local orthonormal frame 
(see Section 3.1) then 


UsUs 


^ UikUik 


which is the desired result, see also for instance [H formula (3.6)]. The Hessian is then replaced by the 
expression in (14.71) and inserted in (14.61) . This leads to the inequality 

(4.8) Ai f |Vu|^d/r<- f — iVul^dcr— f Ric (Vu, Vu) d/r-I-a f |VupdtT. 

Jn 2 Jqq ov Jq Jqq 

Now by (13.311) and (I3.35|) . taking into account that = —au on dVl, we have 

- —iVitP = Il(Vit, Vu) — alVuP -I- auAu + aufiu) + a^(m — 1 )vJH 
2 ov ' ' 

on 9H. Integrating over dVl and substituting into (14.81) we deduce 

(4.9) Ai f \Vuf dp < [ ("iI(Vm, Vu) — q;| Vu|^-b cimAu-I-Q fM/(u)-b a^(m — l)u^id'j dcr 

Jn Jan ^ ' 

— / Ric (Vu, Vu) d/r -b a / |VM|^dcr. 

Jn Jan 

Since dfl is a manifold without boundary we have by the divergence theorem 


/ uAuda = — / IVu^dtr. 

/ an Jan 


On 90 there holds 


|Vu|2 = |Vup -b 


du 


dv 


= |Vu|" + aV. 


Substitution into (14.91) leads to 

Ai f |Vu|^d^ < f ("ll(Vu, Vu)-b a|Vu|^ — 2Qf|VM|^-b a!u/(M)-b a^(m — l)u^7d') dcr 

Jn Jan ^ 2 

— / Ric (Vm, Vm) d/j, 

Jn 

= J (ll(\7u,Vu) — a\Vu\'^ + a^u'^ + auf{u) + a^{m — da 

— / Ric (Vm, Vm) d/i 
Jn 


which completes the proof. 


□ 


We are now in position to state our main result. 
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Theorem 4.4. Let u G C^(fl) be a solution of (|4.1D with f G . Assume that Ric > 0 in fl and that 
for every p G dQ the quadratic form II — ag on the tangent space Tp(dQ), where 5 is the restriction of 
the metric g on Tp(dLl), is nonpositive definite. If in addition 


(4.10) 

then u is unstable. 


/ + auf{u) + a^{m — l)u^H da < 0, 

Jdn 


Proof of Theorem Using (14.51) it is immediate to see that under our assumptions Ai as defined in 
formula (14.21) is strictly negative, so that u is an unstable solution of (14.1|) . □ 

Next we extend this result to the case where condition (14.101) is relaxed relaxed. 


Theorem 4.5. Assume that all assumptions of Theorem 
replaced by 


hold except for condition (I4.10p which is 


(4.11) 


/ + auf{u) + a^(rn — l)u^H da < 0. 

Jan 


Suppose moreover that u ^ 0 and 

(i) a > 0 , or 

(ii) a < 0 and u does not change sign, i.e. either u >0 or u < 0 on LI. 

Then u is unstable. 

Proof of Theorem EH We want to show that under our assumptions Ai as defined in formula (14.21) is 
strictly negative, so that u is an unstable solution of dm. 

We hrst note that u cannot be constant on Ll because the only constant solution of dm is u = 0 . 
Thus, since |Vm| ^ 0 in U it follows immediately from (14.51) and our assumptions that Ai < 0. Let a > 0 
and suppose that Ai = 0. Then |Vm| is a minimizer of the Rayleigh quotient given in (14.2|) . Hence |Vm| 
is a nontrivial eigenfunction associated to the eigenvalue Ai = 0. By the strong maximum principle we 
must have |Vm| > 0 in LI, so that u does not have any critical point in LI. Since LI is compact, u must 
achieve its absolute maximum over H at a point p G dLl and its absolute minimum at a point q G dLl. 
By the Robin boundary conditions and since a > 0 we have 


, , I du, , 
u(,P) = --^(p) < 0, 


, , 1 du , , 

w(<?) =-^(<?) > 0, 


a dv ’ a dv 

so that for every x G LI we have 

0 < u{q) < u(x) < u{p) < 0 , 

which contradicts our assumption u ^ 0. Then Ai < 0, and hence u in unstable. 

Assume now a < 0 and that m > 0 in U. If we assume by contradiction that Ai = 0, arguing as above 
we see that u must achieve its absolute minimum over U at a point q G dLl, where there holds 


Hence we see that 


0 < u{q) = --^{q) < 0 . 
a ou 


<(i) = = 


Since q G dLl is a minimum point for u, all tangential derivatives of u must vanish at q, so that '\/u{q) = 0. 
We conclude that jVuKg) = 0, and hence q is an absolute minimum point for |Vu|. Since |Vm| is an 
eigenfunction of problem dsn associated to the eigenvalue Ai = 0 and since by the strong maximum 
principle we have |Vu| > 0 in Ll, we conclude by the Hopf lemma that 

d 


(4.12) 


dv 


|Vw|(g) <0. 


On the other hand, by the Robin boundary condition in (14.41) . we must have 

■^\L/u\{q) = -a|VM|(g) = 0, 
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which contradicts (14.121) . Thus we have Ai < 0 and u is unstable. 

The case that a < 0 and u < 0 in can be treated in similar way and the proof will thus be 
omitted. □ 

Remark 4.6. Note that the condition II — ag < 0 immediately implies that H = Tr(II) < a on 
dri. Thus, under the above assumptions, condition (14.111) is automatically satisfied if 

1) a > 0 and tf{t) < for every t S K, or 

2) a < 0 and tf{t) > for every t S K. 

Remark 4.7. Notice that if Vm = 0 on dfl, i.e. if u is constant on 
then the hypothesis 

II — ag <0 in Tp{d^l) for every p S 
can be dropped in both Theorems 14.41 and 14.51 . 

Example 4.8. If M is a two-dimensional manifold as in Example 13.31 then the conditions of Theorem 
14.41 (see also Theorem I4.5|l become in view of the computation in Examples 18.31 and 18.51 


each connected component of dfl, 

on 


(4.13) 

Alogp < 0 in 0, 

(4.14) 

Kg -b a > 0 on 90, 

(4.15) 

a — Kn + < 0 on 90. 

au 


Observe that the conditions (14.141) and (I4.15|) coincide with (Cl) and (C2) in Section 2. Condition 
(14.181) is satisfied for a sphere. 


Remark 4.9. If Ric > 0 in O and II — ag < 0 on dft we obtain estimates for stable solutions. In this 
case Ai > 0. Consequently 


fiu) 


> — [a -I- (m — 


Consider the following two examples. 


1 )H] somewhere on 90. 


1. Let f{u) = Xiu, so that u is a solution of Au -I- Aim = 0 in O, -^u + au = 0 on 90 with a > 0. 
Then 

Ai >-[a'^ + {m- I)i7inaxa]- 

2. Let /(m) = —c^M -I- \u\P~^u, so that m is a solution of Am — c^m -b |m|P“^m = 0 in O, ^m -I- om = 0 
on 90 with p > 1, a > 0. Then 


max |m|^ ^ > —\a^ -I- (m — l)i7maxa] + (?■ 
n 

We conclude the section with a Barta type inequality that gives a sufficient condition for stability and 
which will be used in Section [SJ 

Lemma 4.10. Let v be a stationary solution of problem SH) in O. Let there exist a function w G 
C^(0) n (7^(0) such that w > 0 in LI and 

Aw + f'{v)w <0 in Ll 

< 

dw 

— -h aw > 0 

Ol' 


Then v is asymptotically stable. 


on 90. 
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Proof. Let Ai be the smallest eigenvalue of (IE3I and let ipi be the corresponding eigenfunction. We have 
0 > J (pi{Aw + f{v)w}dV = J w{Aipi + f{v)(pi} dfi + J 

/ wipidfi + a / {wipi — (fiw) dv = —Ai / wipida. 

J Q. J an j n 


> —Ai 

/n Jan 

Therefore Ai > 0, thus the conclusion follows. 


□ 


5. Surfaces of revolution in 

A surface of revolution Sjjj in is obtained by rotating around the z-axis a simple, regular plane curve 
r —>■ (i/;(r),x(r)) {r € I = [ri,r 2 ]; ri < r 2 ) with > 0 in (ri,r 2 ). Therefore it admits a parametrization 
of the form 


(5.1) 


X = if{r) cos 0 
y = ipir) sin0, 


((r,6>) e [ri,r2] x [0,27r)) 


z = x{r). 

We can always assume that + (xO^ = 1 in /. Moreover, we suppose that > 0, ifir^) > 0, thus 
dS^ = {(V'(r)cos0,'(/)(r)sin6>,x(r)) |r e {ri,r2},0 e [0,27r)} . 

A surface of revolution in R^ (with parameterization (15.11) 1 is a 2-dimensional Riemannian manifold 
with metric 

ds^ = dr"^ +il}^{r)d9^ . 

In the coordinates (r, 0) (r £ (ri,r2),0 £ (0,27r)) the Laplace-Beltrami operator on is expressed as 


(5.2) 


d'^u du 1 d'^u 

dr^ -ip dr 'ip'^ 89'^ 


A direct calculation shows that the Ricci (Gaussian) curvature of is 

r{r) 


(5.3) 


R{r) = -■ 


ip{r) 


(r e {ri^r^)) ■ 


Observe that it does not depend on the direction X, nor on the angle 9. This is in accordance with the 
fact that on 2-dimensional surfaces the Ricci curvature is independent of the direction and coincides with 

if' 

the Gaussian curvature. Let us also point out for further references that the quantity — represents the 


geodesic curvature kg of the parallel circles r = constant on . 


tp 


5.1. Instability. Let 

O := {{tfir) cos 9, ip{r) sin0, x(r)) | (r, 9) £ [0, a] X (0, 27r]} 

be an annular domain on a surface of revolution with parametrization (15.111 (ci < 0 < a < r 2 ). Note 
that dfl is made of the two geodesic circles: 

Co := {(i/’(0)cos6»,'!/)(0)sin6»,x(0)) |6» £ (0,27r]}, 

Ca ■■= {{ip{a)cos9,ip{a)sm9,x{a))\9 € . 

For the sake of simplicity we assume that x^(0) > 0, x'i^) > 0 ■ 

Let us start with a simple observation concerning non radial equilibrium solutions (see also [5], |22) 
for the case a = 0). 
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Proposition 5.1. Every equilibrium solution v of problem USD, which depends on the angle 6, is 
unstable. 


Proof. By (15.21) ?; is a solution of 


(5.4) 


1 

Vrr + -^Vr + + f {v) = 0 in (ri,r 2 ) X (0,27r). 


If we differentiate this equation with respect to 6 we see that ve is an eigenfunction of (ESI) and that 
A = 0 is the corresponding eigenvalue. The function vg changes sign and therefore it cannot be the 
eigenfunction corresponding to the lowest eigenvalue. Hence Ai < 0, which establishes the assertion. 
□ 


From now on let v{r) be a radial stationary solution. If we differentiate dsa with respect to r we 


get, setting ' := 


Av' + = 0- 


Multiplication by v' and integration over H yields 

— f (u")^ dV + f f'{v){v')'^ dV + f v'v" ds 

Jn Jq Jcc 


— V V 

JCo 


/ 

Jn 


ivrdV = 0 . 


Hence 

(5.5) 


Ai J dV < J {v'Y dV + La{v'{a)v"{a) + a[v'{a)f} - Lo{v'{Q)v"{Q) - a[v'{Q)f) , 


where 


Lo := 27n/'(0), La := 2TTtj}{a). 

Note that for p G Co, X G Tp(Co), one has X = 7 ^ (for some 7 G M) and thus (see, e.g., [5D]) 

H(X,X) = - 7 V(a)f/''(a) ■ 


Hence 

(5.6) 

Similarly, for any q G Co one has 

(5.7) 

Thus, also using dsa, 


H = Ha 


tpia) ■ 


H = Ho = 


im 

V’(O) ■ 


(5.8) 

La{v'{a)v''{a) + a[u'(a)]2} - Lo{v'{0)v"[0) - a[u'(0)]2} 

= La [Haa^v{af' + av{a)f{v{a)) + a^v{af] + Lq [Hoa'^v{0)'^ + av{0)f{v{0)) + a^u(O)^] 


Therefore, we have the next result. 


Theorem 5.2. Suppose that v is a radial stationary solution of E3- If 

<-) -(D.. 

and 


fj" 

"0 


> 0 


(0,a), 


(5.10) La [Haa^v{a)^ + av{a)f{v{a)) + a^u(a)^] + Lq [Hoa^v{0)^ + av{0)f{v{0)) + a^u(O)^] < 0 , 
then V is unstable. 


The assumption (15.911 has a geometrical meaning in the sense that 




(tp) ^ ~ ^ [ri,r2]), 

where Kg{r) is the geodesic curvature of the circles r =const. 
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5.2. Existence of stable solutions (patterns). If the variational problem 

min £(v,D,) 

Dewi'2(n) 

is solvable then the minimizer is stable. Hence for positive a and large classes of nonlinearities this is 
often the case. For Neumann and Robin boundary conditions with negative a the minimum does in 
general not exist. 

In this section we construct on surfaces for which condition (Ell) is violated a problem with negative 
a possessing a stable solutions satisfying the boundary condition (15.1011 . 


Theorem 5.3. If, for some R S (0,a), 



(i?) >0, 


then there exists f € C^(IR.),a < 0 such that problem dOJ admits a stationary asymptotically stable 
solution which satisfies (15.101) . 


In the proof we follow the arguments used in [5] for the case 0 = 0 (see also [5^ where a different 
differential operator is treated). Several modifications are needed to adapt those proofs to our problem; 
they are summarized in Remark 15.61 

Take Rq < Ri < R2 < R3 in a neighborhood of R. Since ip G C'^(-f), we can choose Rq and R 3 such 
that 


(5.11) 



(s) > 0 


for any s G [i ?0 7 .R 3 ] ■ 


Let zi 


(5.12) 


where 

(5.13) 


21 ( 5 ) be the solution of the Cauchy problem 

111 

-Bz = 0 


m' 


ip 


[0,Ri 


z(0)=0, z'(0) = l, 


B > B := max 

[O.a] 



Similarly for any /3 > 0, let Z 2 = ^ 2 ( 5 ) be the solution of the Cauchy problem 


(5.14) 


WV 

. V- . 


-Bz = 0 


in (i? 2 , a] 


[ z{a) = f , z'{a) = —l. 

If necessary we shall write zi = ^i(s, B), Z 2 = Z 2 {s, B, f) to stress the dependence of the solution on the 
parameters B and /3. 


Lemma 5.4. The solution zi of problem (15.121) has the following properties: 
(i) zi > 0 m (0, i?i); 

{ii) zi{-,B) is increasing in [0, i?i) for any B > B; 

(in) zi{r,-) is increasing on {B,oo) for any r in (0,i?i); 

(iv) lim Zi{r,B) = 00 for any r G (0, i?i) . 

B^oo 
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Similarly, for the solution Z2 of problem ( 15 . 131 ) the following hold: 

{i') Z2> P in (i?2, a); 

{a') Z2{-,B) is decreasing in (i?2,a) for any B > B; 

(in') Z2(r, •) is increasing on (B,oo) for any r € (i?2,a); 

{iv') lim Z2{r,B) = c» for any r G (i?2,a) ■ 

B—>-oo 


Proof. The statements concerning zi have been shown in [ 2 ]. Let us show those concerning Z2- 
{i') Assume that there exists r G (i?2,a) such that 

Z2(^ = P, ^2(3) > P for any s G (r, a). 

Then for some f G (r, a) we have 

Z 2 {r) = inaxz 2 > P, z'^if) = 0, z'({f) < 0. 

[r,a] 

So, 


if 


(r) - Bz2{r) = z'f{r) + + 




Z 2 {r) < 0 . 


This contradicts the definition of Z2, hence the claim follows. 


{ii') Suppose on the contrary that there exist r G {R2,a) such that 

( 5 . 15 ) Z2{r) < 0 for any r G {r,a), Z2(r) = 0 => z'fir) < 0 . 

On the other hand, we have 


Z2ir) =-jZ2ir)- 


if 


-B 


Z2{r) > 0 


since by (?) Z2{r) > 0 . This is a contradiction, thus Z2 is increasing in (i?2, a). 

{in') Let B2 > Bi > B. Set Ci(?') := -22(1’; ^i), (2(1') := Z2{r]B2). Then w{r) := (i - (2 solves 

= BiCi - B2C2 < B2C1 - B2C2 = B2W in [R2,a) 


w{a) = 0 , w'{a) = 0 . 


Therefore, w satishes 


w" + 


+ -^2|w<0 in [R2,a) 


w{a) = 0 , w'{a) = 0 . 

Hence, it is easily seen that w < 0 in [i?2,a), so 

Z2{r,Bi) < Z2{r,B2) for any r S [i?2,a), 

thus the claim follows. 


{iv') Fix any Bi > B. Integrating the differential equation in (I 5 . 14 |) and using {ii') we get for any 
r G [i?2,a) and B > Bi: 

1 


( 5 . 16 ) 


22(i’,H) = if{T) J Z2{t,B)dtdT + Pif{a) + 1 +j 'ff{T)d- 

“ I + 1 + -^^/? y '0(r)dr| . 


t)- > 


The claim follows by letting i? —>• 00. 


□ 
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Define 


2:1 (r) 

if r e [0,i?i) , 

zsir) 

if r e [Ri , R 2 ], 

Z2{r) 

if r G (i? 2 , a]; 


(5.17) 


here zg is any positive smooth function such that z is smooth at the points r = Ri, r = R 2 . By its 
definition and Lemma 14.101 ft), the function 2 ; is smooth in [0,a] and 

(5.18) z>0 in (0,a), z{0) = 0, z{a) = P. 

Clearly, 2 : depends on the choice of the parameter /3; to highlight this we write z = zp, if it is needed. 


Lemma 5.5. Let ,0 > 0, let the function z = zp be defined by (15.171) . Then there exists f G C^{ 
that the function 


such 


(5.19) 


Z{r) := f z{s)ds (re[0,a]) 
Jo 


is a stationary solution of problem (nm, which satisfies (15.101) . provided 

P 

In z(.r)dr ■ 


Proof. Since 2 : > 0 in (0, a), the function u = Z{r) is increasing in (0,a). Denote by r = Z ^(u) the 
inverse function, then define 


(5.20) 


f{u) := < 


-Bu-l 

d{f}[Z~^ {u)]z[Z~^ {u)]) / du 


if tt < 0 . 


if 0 < u < Z(a), 


_ -Bu + BZ{a) + l-P^ iiu>Z{a). 

In order to guarantee that / S C'^(R) we have to prove that / is smooth at m = 0 and u = Z{a). The 
smoothness at tt = Z (a) will follow, if we can show that 


(5.21) 


/(tt) = —Bu + BZ{a) + 1 — 


.pm 

V’(a) 


for any tt € (2'(i?2), Z{a)] . 


For that purpose, let us integrate the differential equation in (15.141) on (r, a) for any fixed r G (R 2 ,a). 
We obtain 


(5.22) 


i^^(r) = B[Z(r) - Z(a)] + z'(a) + z(a) for any r G (i?2, a). 


On the other hand, it is easily seen that 

(5.23) 

for any r G (0,a). Therefore, by ()5.22ll - (I5.23|) we have 

(5.24) 


f[Zir)] = - —i^)= - 


>'( 0 ) 


f[Z{r)] = —BZ{r) + BZ{a) + I — P for any rG(i? 2 ,a). 

y^[cL) 


Since Z is increasing, (15.241) holds for tt = Z{r) in {Z{R 2 ), Z{a)]. Similarly it is seen that 
(5.25) /(tt) = —Bu — 1 for any tt G [0, Z{Ri )), 

which implies the smoothness of / at tt = 0. Hence / G C'^(]R). 
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Observe that due to (15.181) and our choice of a we have 

Z'{Q) + aZ{0) = z(0) = 0, Z'{a) + aZ{a) = P + a f z{r)dr = 0. 

Jo 

Note that by (I5.23L Z solves the differential equation in (II.21) . Therefore, Z is a stationary solution to 
problem dOD. Moreover, due to (15.811 it is easily seen that Z ^ 0 satisfies (j5.10|) . Then the conclusion 
follows. □ 


Now we are in position to prove Theorem 15.31 

Proof of Theorem Let Z be the stationary non constant solution of problem with the function 
/ defined by (15.201) of Lemma 15.51 Define 


(5.26) 


'j{r) := < 


z(r) - miz(Ro)(r - Riy 

ifr G [0,Ri), 

z(r) 

if r G [i?i,R 2 ], 

z(r) + m 2 z(R 3 )(r - i? 2 )^' 

if r G (i ?2 , a] , 


with constants mi € (0,cx)),m2 G (0,oo),l G N,l odd, that will be chosen later. Observe that w > 0 
in [0,a]. Furthermore, recall that z, and hence w, depend on the parameter B in problems (15.121) and 

(Oil) . 

Without loss of generality, we can suppose > 0 in (0, a) (see (15.11) 1. Therefore, 

9w , , ,, , dw , , ,, , 

_(0) = -u;(0), —ia) = w{a). 

Next we shall prove the following 

Claim: There exist mi > 0, m 2 > 0,1 G N, / odd and B > 0 satisfying (15.131) such that 

{ifw'Y 

(5.27) 




+ f'{Z)w < 0 


in ( 0 , a) 


rc'(O) — aty(O) < 0 , w'{a) + aw{a) > 0 . 

In order to establish the first inequality in (I5.27L we think of the interval (0,a) as the disjoint union 
(0, a) = (0, i?i) U [i?i,i? 2 ] 0 (i? 2 ,a)- Recall that by definition z = zi in (0,i?i) and z = Z 2 in (i? 2 ,a)- 
Observe that for any r G (0, i?i) U (i? 2 , a) we have by (I5.2ip and ()5.24l) 

ffZir)) = -B 

and by (15.121) and (I5.14p 

tell (3)1)7 

V' V V’/ 

This together with the definition of w yields for any r G (0,i?i) 


^^-'UnZ)ir=^^-B. = - 


(5.28) 


if 


if 


if 


(r) 


+ miz{Ro){Ri - r)^' ^ 3/(31 - 1) + 3/ ( — ) (r)(r - Ri) - B{r - i?i)^ 


if' 


Let us prove that the right-hand side of the above expression is negative in (0,i?i) = (0,i?o) U [i?o,7?i)- 
For this purpose observe that: 

• in ( 0 ,i?o) there holds 




(r) < 


if ^ 


(r) < Bz(Ro) 


(with B defined in (I5.13|) L since z = Zi is increasing by Lemma [501 - ("bll: 
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in [i?Oj^i) we have 


where 




(r) < -Bz{r) < -BziRo ), 


R ■ 

^ := mm — . 

[Ro,Ri] J 

Note that S > 0 by (15.111) : moreover, Lemma [53 ]-(m) has been used again. 


By the above remarks, we have in (0, Rq) 


(5.29) - 


if 


where 






(r) + miz(i?o)(i?i - 3^(3/ - 1) + 3/ j (r)(r - Ri) - B{r - R^f 

< z(Ro) {B + mi{Ri - Rof^-^ [3l{CRi + 31-1)- B{Ro - Ri)^] } , 

3/(3/- l + Ci?!) 


< 


B > 


{Ri - Ro)^ ’ 


C := max 

[0,a] 


Similarly, in we have: 


(5.30) - 


ij 


(r) + miz(i?o)(i?i - r) 


31-2 




31(31 - 1) + 31 ( ^ ) ir)(r - Ri) - B(r - i?i)" 


< 


< z{Rq) [-B + 3lmiRf-^ (31-1 + CRi)] . 

It is easily seen that the right-hand sides of inequalities (I5.29|) . (15.3011 are both negative if we further 
require that 

;B-k31mi(i?i-i?o)^'”^(C'i?i-k31-l) B 


^ mi(Ri-Ro)^^ 

Then from (15.281) we obtain that 

(ijjw'y 


and 0 < mi < 


3lRy-‘^(CRi+3l-l) 


(5.31) 


y, 


+ f(Z)w<Q in (0,i?i]. 


It is similarly seen that, for m 2 > 0 small enough and B > 0 sufficiently large, 


(5.32) 




-|- f'{Z)w <0 in (i? 2 ,o)- 


Now consider the interval [i?i, i? 2 ]- Since Z is a stationary solution of problem (11.211 . in [i?i, i? 2 ] there 
holds 

Z" + ^Z' + f(Z) = 0. 

ip 

Deriving the above equality and recalling that Z' = z we obtain 

z"+'^z'+ f'{Z)z = - z in [i?i,i? 2 ]. 

The right-hand side of the above equality is negative in [i?i,i? 2 ] by inequality (I5.11|l . thus 
(5.33) 




iP 


iP 


From (I5.32I1 - (I5.33I1 we conclude that the first inequality of (15.2711 is satisfied. It remains to prove the 
inequalities —w'(0) + q:u;(0) > 0, w'{a) -I- aw{a) > 0. For this purpose, note that in view of (15.1611 we 
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can infer that there exist two constants Co > 0, B* > B such that 

f z{r\B) dr > CqB for any B > B* . 

Jo 

Thus, 

(5.34) a = — ra > - 7 ^ foranyS>B*. 

Jp z[r)dr Oq-d 

Hence, in view of (15.341) . {iv) and {iv'), choosing B > B* large enough and I > max{i?i, a — R 2 }, 
we obtain 

—m;'(0) + aw{0) = —1 + miz{Ro)Ri’‘~^{31 + aRj) > 0 , 

and 

w'{a) + aw{a) = —1 + a/3 + m 2 z{Ro){a — i? 2 )^*~^[ 3 Z + (a — i? 2 )a] > 0 . 

This completes the proof of the Claim. Observe that (I5.10|) is satisfied. Then by Lemmas 15.51 and 
14.101 the function Z is a stable stationary solution of problem (11.21) with / given by (I5.20|) . Then the 
conclusion follows. □ 

Remark 5.6. Note that the construction of z and / are similar to that in [5]. However, in [5] we had 
/3 = 0; instead now we need /3 > 0. Moreover, in the proof of the result in [2] analogous to Theorem 15.31 
we had Z = 1 in ()5.26p . 


6 . Further examples 

6.1. Spherically symmetric manifolds. We start recalling some basic notions on spherically sym¬ 
metric manifolds. Let M be a complete Riemannian manifold. Let us fix a point o £ M and denote by 
Cut(o) the cut locus of o. For any x £ M \ [Cut(o) U {o}], one can define the polar coordinates with 
respect to o, see e.g. [5]. Namely, for any point x £ M \ [Cut(o) U {o}] there correspond a polar radius 
r{x) := dist{x, o) and a polar angle 9 £ such that the shortest geodesic from o to x starts at o with 

the direction 9 in the tangent space TqM. Since we can identify TqM with M™, 9 can be regarded as a 
point of 

The Riemannian metric in M \ [Cut(o) U {o}] in polar coordinates reads 

ds"^ = dr"^ + A,j{r,9)d9^d9^, 

where {9^,... ,0™“^) are coordinates in and (Hy) is a positive definite matrix. It is not difficult 

to see that the Laplace-Beltrami operator in polar coordinates has the form 

(6.1) A = ^+.F(r,e)|; + As„ 

where R{r,9) := ^(log Y^H(r, d)), A{r,9) := det{Aij{r,9)), As^ is the Laplace-Beltrami operator on 
the submanifold Sr ■= dB{o,r) \ Cut(o). 

M is a manifold with a pole, if it has a point o £ M with Cut(o) = 0. The point o is called pole and 
the polar coordinates (r, 9) are defined in M \ {o}. 

A manifold with a pole is a spherically symmetric manifold or a model, if the Riemannian metric is 
given by 

( 6 . 2 ) ds^ = dr^ + (j)^{r)d9^, 

where d9^ is the standard metric in 8™“^, and 

(6.3) e A := {/ e C°"((0, 00 )) n Ci([0, CX))) : /'(O) = 1, /(O) = 0, / > 0 in (0, 00 )}. 

In this case, we write M = furthermore, we have yjA{r, 9) = (fC~^{r), so the boundary area of the 
geodesic sphere dSn is computed by 


S{R) = u:m4>^-\R) 
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uJm being the area of the unit sphere in W^. Also, the volume of the ball is given by 

Vol(Bi^(o))= 

Jo 


Moreover we have 


A “ ^ ^ 1a 


A 52 S' d 1 , 

A — „ „ H— 1 rrrAsm-i 


or equivalently 

^ ^ dr'^ ' S dr ' (j>^ 

where Agm-i is the Laplace-Beltrami operator in Note that similarly to (15.61) and (15.711 one can 

compute the mean curvature of dBp(p) in the radial direction ^ as follows 


(6.5) 


d)' (r) 

H(r) := ——— for each r > 0. 
d>{r) 


Observe that for (j){r) = r, M = M™, for (j){r) = sinhr, M is the m—dimensional hyperbolic space 
H™, while for (f){r) = sinr (r S [0,7r)) we have the m—dimensional sphere S’" C (see [5]) . 

For any x £ M \ [Cut(o) U {o}], denote by RiCo(x) the Ricci curvature at x in the direction If 
M = is a model manifold, then for any x = {r,6) £ M \ { 0 } 


( 6 . 6 ) 


RiCo(a;) = —(m — !)■ 


(j){r) 


Now we discuss the stability of radial solutions of problem (11.211 with O := Bn^o) \ Br{o) C for 
each 0 < r < R. 

In view of (|6.4D and (15.4L setting S = tp, the same results as in Section [5] hold. Indeed, we have the 
following theorem. 


Theorem 6.1. Let 17 := Bb.{o) \ Bp{o) C with 0 < p < R. 


(i) Suppose that v is a radial stationary solution of If 


(6.7) 


- ^ =-^ + bT >0 (P’R) 


and 

d,{R){v'{R)v"iR) + a[v\R)n - 0(p)K(p)u"(p) - a[v'ip)n 

( 6 . 8 ) = (piR) [a^(m — l)iJ(i?)u(i?)^ + av{R)f{v{R)) + a^v{R)‘^] 

+ Hp) [a^im - l)H{p)v{pf + avip)f{v{p)) + a^v{pf] < 0 , 

then V is unstable. 

(ii) If for some R £ {p, R) 

( 6 . 9 ) (i?)> 0 , 

then there exists f £ (K), a < 0 such that problem (O admits a stationary asymptotically 

stable solution which satisfies (IS3). 


Note that, in view of (1531) and (16.6L the inequalities (16.71) and (16.91) have a geometrical meaning. 
Indeed, (EZl) is equivalent to the following requirement 

RiCo(x) > —(m — l)[H{r)]^ for any x = {r,0) £ LI, 


and similarly for (ESI). 
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6 . 2 . 


Straight cylinder in 


A straight cylinder C in is parameterized as follows 

X = 'ljj{t) 


( 6 . 10 ) 


< y=x{t) ((i,s) G X [si,S2]) 


z =s, 

where t i->- {ip{t),x{t),0) is a simple, regular, closed plane curve {t G [ti,t 2 ]; < ^ 2 )- We suppose that 

[4’'{t)]'^ + [x!{i)Y = 1 ^ G [ti, < 2 ] • It is easily seen that, for all p G C, A G TpC, 

Ric(A,A) = 0; 


furthermore, since the second fundamental form of dC with respect to the embedding dC ^ C is identi¬ 
cally zero, we also have that its mean curvature identically vanishes. 

Note that 

Au{t,s) = Utt(t,s) +Uss(t,s) . 

Then, by a similar argument to that of Proposition 15.11 one can see that any stable solution of problem 
(EH) must depend only on the variable s. 

Now, consider a solution u = u(s) of problem Thus, using the same notation as in Section |4l 

we have Vit = -^u = 0. Hence, from the same arguments used in the proofs of Proposition 14.41 and of 
Theorem 14.51 we can infer that if 

i) / [a^u^ -I- auf{u)] dcr < 0 , or 

JdC 

ii) a > 0 and / [a^u^ + auf{u)] da < 0 , or 

JdC 

iii) a < 0, [a^u^ + auf{u)] da < 0 and u does not change sign, 

JdC 

then u is not stable (see also Remark 14.71) . 
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